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ABSTRACT 

Methane is the most abundant simple organic molecule in the outer solar system 
bodies. In addition to being a gaseous constituent of the atmospheres of the Jovian 
planets and Titan, it is present in the solid form as a constituent of icy surfaces 
such as those of Triton and Pluto, and as cloud condensate in the atmospheres of 
Titan, Uranus, and Neptune. It is expected in the liquid form as a constituent of the 
ocean of Titan. Cometary ices also contain solid methane. The optical constants for 
both solid and liquid phases of CH^ for a wide temperature range are needed for 
radiative transfer calculations, for studies of reflection from surfaces, and for 
modeling of emission in the far infrared and microwave regions. The astronomically 
important visual to near infrared measurements of solid methane optical constants 
are conspicuously absent from the literature. 

We present preliminary results on the optical constants of solid methane for 
the 0.4 nm to 2.6 urn region. Deposition onto a substrate at 10°K produces glassy 
(semi-amorphous) material. Annealing this material at ~33°K for ~1 hour results in a 
crystalline material as seen by sharper, more structured bands and negligible 
background extinction due to scattering. We report A: for both the amorphous and the 
crystalline (annealed) states. Typical values (at absorption maxima) are in the 10' 3 
to 10 range. Below A = 1.1 /tm the bands are too weak to be detected by transmission 
through the films < 215 /im in thickness, employed in our studies to date. 

Using our previously measured values of the real part of the refractive index, 
n, of liquid methane at 110°K [Bull. Am. Phys. Soc. 31, 700 (1986)] we compute n for 
solid methane using the Lorentz-Lorenz relationship. Values vary monotonically from 
n = 1.304 at A = 0.4 n m to n = 1.296 at A = 2.0 fim. An independent, direct determination 
using interference fringes at two angles of incidence finds n = 1.31 ± 0.02 for the 
solid film at A = 0.4-0. 7 /mi, a value confirmed by laser interferometry. 

Work is in progress to extend the measurements of optical constants n and k for 
liquid and solid to both shorter and longer wavelengths, eventually providing a 
complete optical constants database for condensed CH 4 . 
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INTRODUCTION 


The absorption, emission, reflection and scattering of radiation by at- 
mospheres and surfaces depends on the vertical and horizontal distribution of 
material, on particle size and shape, and most fundamentally on the optical charac- 
teristics of the materials present. These optical characteristics can be represent- 
ed by the complex refractive index. The real and imaginary parts are commonly called 
the optical constants n and k , and can be determined in the laboratory by a variety of 
techniques. 1 CH 4 optical constants are required in radiative transfer calculations 
modeling light scattering and the thermal structure of planetary atmospheres, and in 
studies of reflection from particulate surfaces containing CH 4 ice and frost, such as 
those of Triton and Pluto; and also in the calculation of the energy deposition 
profile and penetration depth for charged particles impinging on icy outer solar 
system surfaces. These quantities can be calculated from the "stopping power": 



where e = e 1 + e 2 , = n 2 - k z , e 2 = 2 nk, n and k are the real and imaginary part of the 

refractive index, amd u> is the frequency. 

We have been engaged for the last several years in the laboratory determination 
of complex refractive indices of materials of planetary and astronomical interest. 

The optical constants of Titan tholin 1 produced from ~ 90% N 2 and 10% CH 4 in a plasma 
discharge similar to that produced by Saturnian magnetospheric charged particles 
which are intercepted by Titan’s upper atmosphere, were utilized in radiation 
transfer calculations, 2 ' 3 ' 4 producing, for example, reflection spectra 5 and thermal 
structure 4 of Titan’s atmosphere consistent with the ground based spectroscopic 
observations, 6 with Pioneer 11 polarimetric observations; 7 and with Voyager photome- 
tric observations. 8 ' 9 Successful applications were also made with optical constants 
of H 2 /He/CH 4 tholin appropriate to Neptune and the Uranian troposphere. 10 One of our 
current emphases is the determination of optical properties of condensed hydrocar- 
bons and the effect of irradiation on these properties. Simple hydrocarbons such as 
methane (CH 4 ), ethane (C 2 H 6 ), and acetylene (C 2 H 2 ) are prime candidates for such 
studies because of their importance on several planets and their icy satellites. For 
example, 

(a) Titan’s surface probably has lakes or oceans of CH 4 and C 2 H 6 , with several 
other minor species. 11 ' 12 ' 13 

(b) C 2 H 6 , and C 2 H 2 form condensate hazes in the stratospheres of Uranus and 
Neptune. 14 

(c) CH 4 clouds are present on Uranus and Neptune. 15 


328 



(d) Condensed CH 4 is present on both Triton 16 and Pluto 17,18 and 

(e) The presence of CH 4 in cometary comae 19 suggests the presence of condensed 
CH 4 and probably other hydrocarbon ices there. 

A critical survey of the literature reveals that the availability of optical 
constants data for these materials in the visible-near IR region ranges from sparse 
to totally absent. We present below some preliminary results of the measurements of 
optical constants of solid CH 4 in the 1.1 /xm to 2.6 nm range, and calculations of k and 
a tabulation of n throughout the visible and near-infrared. 

EXPERIMENTAL DETAILS AND RESULTS 
Sample Preparation 

The solid films of methane were made by pulse deposition of gaseous methane on 
a substrate attached to the cold finger of a commercial cryocooler 20 as shown in 
Figure 1 . An externally controlled heater wound around the cold finger of the 
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Fig. 1. Experimental set-up to measure transmission through 
methane film at 10°K 

cryocooler enabled us to control the sample temperature to within ± 0.5°K. Both KRS-5 
(for transmission) as well as polished silicon (for reflectance) were employed as 
substrates. Typically, the sample was deposited at 10°K which resulted in amorphous 
films, as evidenced by broad structureless bands in the infrared region (Figure 2). 
Annealing the sample at ~33°K for approximately an hour resulted in crystallinity, as 
evidenced by sharper structured bands and negligible background extinction due to 
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scattering (Figure 3). Several films from 5 /xm to 215 /x m thickness were investigated. 
Spectral Measurements 

The spectra were recorded on a Cary- 14 UV-VIS-NIR spectrometer from 0.35 /xm to 
2.5 /xm with a spectral resolution ~6 cm' 1 . The region 2.0 /xm to 2.6 /xm was also 
measured with a Perkin-Elmer 1800 FTIR spectrometer at ~1 cm’ 1 resolution. 

Data Analysis 

Imaginary part of the refractive index k 

The determination of the imaginary part of the refractive index, k, is deter- 
mined directly by measuring the transmission through a film of thickness t and using 
the relation. 


k = — In 

bnt 


t s a-R ) 
T (1 -R 0 ) 2 


( 1 ) 


where T = transmission, T g = Transmission of the substrate alone. 



and n g = the real part of the refractive index of the substrate, known from indepen- 
dent measurements or from the literature. 21 R = the total reflectance of the film and 
substrate, which can be approximately determined by constructing a baseline through 
points where the film has no absorption bands. Further details are given by Khare et 
al. 1 

The only unknown in Equation (1) is the thickness of the film t which we 
determine from the interference fringes at 0.35 /xm to 0.75 /xm for two angles of 
incidence. 22 Since the absorption in this region is very small, the real part of the 
refractive index does not undergo any appreciable dispersion. Assuming n to be 
constant, the two equations 


It ( n 2 

and 

2t (n 2 


sin 2 ©! ) 1/2 
sin 2 0 2 ) 1/2 


/ ^ 1^2 
/ ^ 3^4 


) 

) 


( 2 ) 


where AN is the number of fringes in the range (A., - A 2 ) and (A 3 - A^) for the two angles 
of incidence 0, and © 2 respectively, can be solved for n and t. The thickness of the 
film can be obtained within about 2-5% accuracy for 5-30 /tm thick films. For the 


331 



thickest film (215 /im) employed in our work, the fringes are too close to be separated 
accurately. For this sample, the thickness was estimated from the integrated band 
intensities of the stronger peaks and from the knowledge of the extinction coeffi- 
cient of these bands measured for thinner films. 

Substituting for t in Eq. (1) yields the k values for crystalline CH 4 tabulated 
in Table 1 and plotted in Figure 3. 

The k values calculated for the amorphous film of methane are plotted in Figure 
2. These results are of much poorer quality than those for crystalline CH 4 because of 
large background scattering. If the background contribution is subtracted from the 
observed data, then the k values for amorphous methane are approximately the same 
magnitude as the corresponding k values for crystalline methane. As we go towards 
wavelengths shortward of A = 1.6 /im (Figure 2) into the visible, the loss of photons 
due to scattering dominates over the loss due to absorption at the peaks to such an 
extent that one finds complete loss of absorption features. The k values at peaks for 
solid glassy amorphous film are higher than those of the corresponding peak values of 
the crystalline film at ~33°K because of the background scattering by the non- 
crystalline sample. 

For crystalline CH 4 (a 215 /im thick sample), scattering is much less of a 
problem, allowing determination of values of k as small as 10' 6 down to a wavelength 
of 1.1 /tm. 

Extinction in solid CH 4 in the VIS-NIR region (0.4- 1.1 /im) is expected to be 
extremely small ( k < 10' 6 ). In Figure 4 we plot k for liquid CH 4 23,24 (solid line) 
measured by Ramaprasad 23 et al. and Patel 24 et al. and the k values computed 25 from 
room-temperature gas data obtained by Giver 26 if taken to a liquid density of 0.45 g 
cm" 3 (dashed line). At low resolution, both the experimentally measured liquid CH 4 k 
values and those computed from gas data, are similar (Figure 4), and therefore 
temperature and density effects for CH 4 are minor in the VIS-NIR (A < 1.1 /im). Since 
CH 4 retains rotational freedom in the solid state, only minor liquid-solid spectral 
differences are expected. 25 Figure 4 therefore provides a good guide for attempts to 
measure k in the solid at visual wavelengths. We are attempting to make optical 
quality solid methane discs varying in thickness from 1 to 10 cm, but Figure 4 would 
indicate that much thicker samples of the order of a meter or so are required to 
obtain k data in this region. This work is continuing at present, along with inves- 
tigations of weak interband absorption and some work on other hydrocarbons. 

The measurement accuracy for the crystalline CH 4 results of Figure 3 and 
Table I is about ± 5%, and is indicated by error bars on the graph. The largest single 
source of error is uncertainty in the correction of absorbance measurements on the 
Cary- 14 spectrometer, which deviate from linearity at the higher end of the 0-1 
absorbance scale. The value of k determined from the data taken on the P.E. 1800 FTIR 
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Imaginary Part of the Refractive Index (k) at Various Wavelengths (A) of a Crystalline Film 

of CH, at ~33°KL 
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are good to around ± 10%, with most of the error coming from the estimate of the 
215 /im-thick sample. 



Wavelength, jj,m 


Fig. 4. Comparison 25 of k of liquid CH^ 23 * 3 ^ (dashed line) with computed k values from 
room temperature gas data'* if taken to a liquid density of 0.45 g cm' 3 (solid line). 
Values are similar everywhere, and nearly indistinguishable below 0.55 fim. 

Real part of the refractive index n 

Using the technique of ellipsometric polarmetry, we have measured the refrac- 
tive indices of liquid methane 27 at ~1 10 o K from A = 0.4 to A = 2.0 yum. Using the 
Lorentz-Lorenz relationship 




n 2 + 2 / n 2 + 2 


= d J d l 


( 3 ) 
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where subscripts / and 5 represent the liquid and solid phases of CH 4 respectively, we 
could estimate the real part of the refractive index knowing the densities of liquid 
and solid CH 4 . The density 28 of liquid CH 4 at 110°K is 0.42441 g cm' 3 . The density of 
solid CH 4 , d s = 0.449 g cm' 3 at the temperature of our solid film, ~20°K. This density 
is obtained by extrapolating the values of Constantino and Daniels. 29 n values calcu- 
lated by using Lorentz- Lorenz relationship [Eq. (3)] are given in Table II and are in 
good agreement with the value of n = 1.31 ± 0.02 calculated from interference fringes 
we obtained during reflection measurements for two angles of incidence (Eq. 2). 


Table II 


Real Part of the Refractive Index n s of Solid Methane Calculated from 
Lorentz-Lorenz Relationship 


A (/zm) 

0.4 

0.6 

0.8 

1.0 

1.2 

2.0 


1.286 

1.283 

1.281 

1.280 

1.280 

1.278 


1.304 

1.301 

1.299 

1.298 

1.298 

1.296 


n for solid methane at A = 0.6328 moi was further confirmed independently by a 
dual angle laser interference technique 30 employed to monitor film thickness and the 
refractive index. Two He-Ne laser beams were specularly reflected off a polished 
silicon crystal which was used as a substrate (Fig. 5). 
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As the methane gas is condensed at a constant rate, two interference patterns 
of different periods are recorded in the reflected laser light as a function of film 
thickness (Fig. 6). If the ratio of separation between adjacent maxima in two 



beams is termed /?, then 


n = ( sin 2 e b - /? 2 sin 2 8 a ) 1/2 (4) 

( 1 - 0 2 ) 1 ' 2 

where © a and © b are 32.6° and 69.8°, the two laser beam incidence angles. The index of 
refraction from the data in Fig. 6 gave n = 1.302 ± 0.002 at A = 0.6328 /mi, confirming 
again that the value obtained by Lorentz-Lorenz relationship [Eq. (3)] is valid for 

ch 4 . 

DISCUSSION AND CONCLUSIONS 

Buie et al. 17 have obtained photometric data of Pluto-Charon mutual events at 
wavelengths 1.53, 1.75, 2.00 and 2.35 /im while Sawyer 31 has obtained data at wave- 
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lengths from 0.5 to 1.0 pm. This data allows the extraction of spectra of Pluto 
alone. Our optical constants from 1.1 to 2.63 pm (Table I and Fig. 3) and future 
determinations down to 0.4 pm are specifically useful in scattering models for Pluto, 
which are needed to determine how much CH 4 absorption arises f rom the surface and how 
much from the atmosphere. 17, 18,31 In a similar way, they can be used to model the 
texture and coverage of CH 4 on Triton. 32 

Fink and Sill 33 have studied the reflection spectra of CH 4 frost at 60°K at a 
resolution of 11 cm' 1 from 1.12 pm to 3.85 pm. They also determined transmission 
spectra through a film 7.1 pm thick, in the mid-IR from 2.63 pm to 20 pm, at 57°K and 
72°K. Optical constants of CH 4 for this longer wavelength range could be estimated 
from that work, but are not computed explicitly here. Our k values are currently 
determined near band maxima for the near-IR 1.1 pm to 2.6 pm region. Attempts to 
measure k for CH 4 at shorter wavelengths (0.4- 1.1 pm) by direct transmission are 
difficult due to scattering from surface imperfections. These surface scattering 
losses dominate the true extinction losses. To overcome this problem, we are 
investigating the possibility of using a double beam configuration, with a short 
pathlength through solid CH 4 in the reference beam and a long pathlength through 
solid CH 4 in the sample beam. This should largely cancel the scattering from the 
surface and thus allow the extraction of true transmission information. 

Since the optical constants of CH 4 have application in a variety of planetary 
and cometary surface and atmospheric environments, we intend to measure optical 
constants for a wide variety of temperatures. Some work with C 2 H 6 and C 2 H 2 has been 
done, and studies of changes in the optical properties of CH 4 and other hydrocarbon 
ices after ultraviolet and charged particle irradiation are also planned. 
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